Stochastic Ballistic Annihilation and Coalescence 
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We study a class of stochastic ballistic annihilation and coalescence models with a binary velocity 
distribution in one dimension. We obtain an exact solution for the density which reveals a universal 
phase diagram for the asymptotic density decay. By universal we mean that all models in the class 
are described by a single phase diagram spanned by two reduced parameters. The phase diagram 
reveals four regimes, two of which contain the previously studied cases of ballistic annihilation. The 
two new phases are a direct consequence of the stochasticity. The solution is obtained through a 
matrix product approach and builds on properties of a (/-deformed harmonic oscillator algebra. 
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Systems of reacting particles are used to model a whole 
gamut of phenomena relevant to fields ranging from 
chemical physics through statistical physics to mathe- 
matical biology. In some applications the particles repre- 
sent chemical or biological species jUI); m other cases 
they are to be interpreted as composite objects such 
as aggregating traffic jams ||. Excitations can also be 
treated as interacting particles, one example being laser- 
induced excitons in certain crystals H. Furthermore, 
domain walls occurring in a number of different contexts 
such as growth and coarsening processes (||| have dy- 
namics with a natural particle interpretation. 

Generally these systems are defined through nonequi- 
librium dynamics. Given such a wide variety of nonequi- 
librium reaction systems, it is natural to ask if they can 
be divided into distinct groups akin to the universality 
classes known for equilibrium systems. 

Two reactions that have been extensively studied are 
single species annihilation (A+A— >0) and coalescence 
(A+A-^A). A particularly striking result is that if the 
reactant motion is diffusive, the two processes belong to 
the same universality class H and the density decay is 
independent of the reaction rate in two dimensions and 
below. Moreover, these diffusive systems have also served 
as prototypes for the development of a variety of theo- 
retical tools ranging from field theoretic renormalization 
group (RG) H , to exact methods in low dimensions j| . 

On the other hand, much less is known about the same 
reactions when the motion is ballistic (deterministic) de- 
spite the relevance of such motion to the modeling of 
growth and coarsening processes |||6| . A seminal model 
was introduced and solved by Elskens and Frisch [jl0| and 
describes pairwise annihilation of oppositely moving par- 
ticles in one dimension (ID). That study was restricted 
to particles that react upon contact with probability one. 
More recently some results have been obtained for sys- 
tems in which the reaction probability is less than one, 
thus introducing stochasticity into the evolution Jll|,|l2| ■ 

In this work we introduce a class of ID stochastic bal- 
listic reaction systems. The class includes both ballistic 



annihilation and coalescence and incorporates as special 
limits the models of |l(J and Jl2] . We obtain an exact so- 
lution for the density decay which reveals a single phase 
diagram common to all combinations of ballistic annihila- 
tion and coalescence. This demonstrates a universality of 
the two processes. The universality is stronger than that 
usually discussed in an RG context and can be likened to 
a law of corresponding states. The phase diagram generi- 
cally comprises four decay regimes in contrast to the two 
previously known JuH . The new phases are a result of 
the stochasticity of the reactions. 

Our exact solution is based on the invariance of certain 
properties of our class of models under change of the 
initial spacing of the particles. As a consequence the long 
time density may be determined exactly through a matrix 
product approach of the type introduced in |l3| ] . We use 
this property, and employ recent results on g-deformed 
algebras to analyse the asymptotic density decay. 

We now define the class of models to be considered. At 
time t = reactants are placed on a line with nearest- 
neighbor distances chosen independently from a contin- 
uous exponential distribution. The unit of length is cho- 
sen so that the initial density is g — 1. Although we con- 
sider here a Poisson initial condition, our methods are ex- 
tendible to more general initial distributions (see below). 
Each particle is assigned a velocity +c (right-moving) or 
— c (left-moving) with probability /r and ft, = 1 — /i? re- 
spectively. Particles move ballistically until two collide, 
at which point one of four outcomes follows, see Fig. [l]: 
the particles pass through each other with probability q; 
the particles coalesce into a left (right) moving particle 
with probability pr/^ (prjji); the particles annihilate with 
probability p(l — t]l — Vr)- Here p = 1 — q is the proba- 
bility that some reaction occurs. 

Before describing our method of solution, we present 
our main results summarized in the phase diagram Fig. |^. 
Two important quantities that emerge are the reduced 
densities /£ = f L (l - m) and f* R = f R {l - r}£). The 
phase diagram is spanned by q, the probability of not 
reacting, and x = SrIIli the ratio of reduced densities. 
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FIG. 1. The possible reactions and their probabilities. 
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FIG. 2. The phase diagram of the model for x — ^- The 
behavior for \ > 1 can be deduced from that with \ < 1 via 
a left-right symmetry (see text). 



For simplicity we consider \ < 1 which results in eventual 
extinction of right-moving particles; \ — 1 is a special 
case where both species die out. The case x > 1 can be 
treated using the symmetry of the model under left-right 
interchange f R <-> f L and vr^ Vl i-e. x -> 1/x- 

The different phases — of which two are regions and two 
are lines in the phase diagram — correspond to four qual- 
itative long-time density decays. When x < q 2 the decay 
is purely exponential of the form g ~ g^ + b\ exp(— ait); 
when x — Q 2 the decay is exponential, multiplied by 
a t~ x / 2 power law Q — Qoo + b 2 exp(— a2 i )/t 1 / 2 ] when 
q 2 <X<l the decay exponential, multiplied by a t~ 3 / 2 
power law g — g^ + 63 exp(— a 3 t)/t 3 / 2 ; finally when 
X = 1 the decay is pure power law g ~ 62 A 1 / 2 - The 
exact expressions for the coefficients a, b and final den- 
sity goo are as given in Table [|. Some intriguing features 
that emerge from the phase diagram Figure ^ are: 

(i) There is universality of ballistic annihilation and 
coalescence. This is manifested by the fact that all the in- 
formation concerning the reactions of a particular model, 
along with the initial densities, are encoded into a single 
parameter x- F° r a generic choice of vr, t)l defining a 
particular annihilation-coalescence model, the same four 
decay regimes are found by varying the initial densities 
or stochasticity parameter q. In this way the universality 
can be considered as a law of corresponding states. 

(n) Two new density decays appear which were not 
anticipated in previous works. The first is the line x = q 2 
and the second the region x < q 2 ■ Thus for a generic 
value of the reaction probability 1 — q, varying the initial 
densities gives rise to four types of asymptotic decay. 

(Hi) The deterministic case q = is non-generic since 
along this line only two of the possible phases are tra- 
versed. For the pure annihilation model (77^ = T)l = 0) 
these phases were found in g(J. Thus we refer to the 



entire region q 2 < x < 1 as the Elskens-Frisch phase. 

(iv) The line x = 1 (equal reduced densities) is non- 
generic since a single, power law, decay regime is found. 
The decay does not depend on the stochasticity q. For 
VR = Vl = this special line corresponds to equal ini- 
tial densities ]l2] j. Our results show that such a special 
line exists for all combinations of annihilation and coales- 
cence. This phase can be understood through the picture 
of jl0| . Density fluctuations in the initial conditions lead 
to trains of left- and right-moving particles: in a length 
~ t the excess particle number is ~ t 1 / 2 which yields 
the i -1 / 2 density decay. At long times, the train size is 
large and so a particle in one train encounters many par- 
ticles in the other and will eventually react making the 
parameter q irrelevant. 

(v) In the two new phases (x < q 2 ) the two species 
decay at unequal rates leaving a non-zero population of 
left-moving particles. This is to be contrasted to the 
Elskens-Frisch phase (x > q 2 ) where the final density of 
one species is non-zero but both species decay at the same 
rate. A simple example of non-equal decays is the case 
Vr, = 0, rjL = 1 (x = 0). Then left-moving particles do 
not decay but simply absorb the right-moving particles 
with probability 1 — q giving g = fne~ 2 ^ 1 ~ q ^ LCt . Our 
results show that, in general, increasing q leads to a non- 
trivial transition at x = 9 2 to a regime where the two 
species have different decay forms. 
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TABLE I. Long time density decays to the asymptotic 
value Qoo — /l(1 — x)- Results for \ > 1 can be obtained via 
the right-left symmetry (see text). K is given by (|Io|). 

We now turn to the method of derivation of the phase 
diagram. The density is given by 



g(t) = fRP ( s H \t) + f L P { t\t) 



(1) 



where Pg L \t) and Pg R \t) are the probabilities that a 
left- and right-moving particle survive up to a time t 
respectively. These two probabilities are related by the 
left-right symmetry noted previously; therefore, if we cal- 
culate P { S R) (t) for all X we can infer Pg L ^ (t). To calculate 

{ R) 

Pa (t) consider the leftmost, right-moving particle in 
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FIG. 3. A configuration and set of trajectories and reac- 
tions for a test particle (shaded and bold line) encountering 
a string of N — 10 particles. Note how changing the spacing 
between, for example, the fifth and sixth particles (indicated 
by dotted lines) alters the time sequence of the reactions but 
not the final survival probability. 

Fig. H which we refer to as a test particle. From the figure 
one can see that the initial spacing of the particles on the 
line does not affect the sequence of possible reactions for 
any given particle, in particular for the test particle (we 
return to this point later). Also note that after a given 
time t, the test particle may only have interacted with 
the N particles initially placed within a distance X = 2ct 
(and to the right) of the chosen particle. These two facts 
imply that the survival probability can be expressed in 
terms of two independent functions. The first is F(N), 
the probability that the test particle survives reactions 
with the N particles initially to its right, and depends 
only on the sequence of the N particles. The second is 
G(N; X), the probability that initially there were exactly 
N particles in a region of size X = let. Explicitly, 



N=0 



(2) 



Thus the problem is reduced to two separate combinato- 
rial problems of calculating G(N; X) and F(N). For the 
Poisson initial conditions G(N;X) = X N e- x /Nl g§. 
In the following we show how the second problem may 
be solved by employing a matrix product approach Jl3| . 

As an example consider a test particle encountering 
the string of reactants depicted in Fig. ||. We claim that 
the probability of the test particle surviving through this 
string may be written as 

(W\RRRLRLLLLR\V) (3) 

where R, L are matrices (or operators) and (W\, \V) vec- 
tors with scalar product (W|V) = 1. Thus we write, in 
order, a matrix R (L) for each right (left) moving particle 
in the initial string. 

We now show that the conditions for an expression 
such as ([?]) to hold for an arbitrary string are 

RL = qLR + p(rj L L + r) R R+[l-ri L -r) R )) (4) 
(W\L=(W\(q+ Pm ) R\V) = \V). (5) 



To understand condition (||), recall that after an inter- 
action between a right-moving and left-moving particle 
there are four possible outcomes (see Fig. 1) correspond- 
ing to the four terms on the right hand side of Eq. ^ with 
probabilities given by the respective coefficients. Using 
(0) , any initial matrix product such as (||) can be reduced 
to a sum of terms of form (W\L S R 1 ^} corresponding to 
all possible final states ensuing from the initial string 
and with coefficients equal to the probabilities of each fi- 
nal state. The test particle will survive such a final state 
and pass through the s left-moving particles with prob- 
ability (q + prjuY . The conditions (^]) ensure that this 
probability is obtained for each possible final state. 

The above approach relies on an important property 
of the system which is invariance of a reaction sequence 
with respect to changes of initial particle spacings. To 
understand this, consider again Fig. |[ By altering the 
initial spacings of the particles, the absolute times at 
which trajectories intersect and reactions may occur (if 
the reactants have survived) may be altered. For exam- 
ple, by increasing the spacing between the fifth and sixth 
particles, the trajectories of the third and fourth parti- 
cles can be made to intersect first. However as we have 
already seen, for any particle, the order of intersections 
it encounters does not change and so the final states and 
probabilities are invariant. This invariance is manifested 
in the matrix product by the fact that the order in which 
we use the reduction rule (Jij) is unimportant i.e. matrix 
multiplication is associative. 

Averaging over all initial strings of length N yields 



F(N) = (W\(f L L + f R R) N \V) 



(6) 



To evaluate F(N) we write R 



Vfyi/fB. a + r) L and 
rjR. One can check from (f||]) that 



a, a* satisfy a q-deformcd harmonic oscillator algebra 



(W\a) = (W\q/Vx 



aa' — qa) a = 1 — q 
a\V) = ^x\V). 



(7) 
(8) 



As is evident from (g) the vectors | W) , | V) are eigenvec- 
tors of a, and are called g-deformed coherent states. The 
explicit form of these eigenvectors is known ]f4| . Using 
the above definitions dol) becomes 



F{N) = {W\[J] R Ji{a + a)) + l-fi-f R \ \V). (9) 

In the deterministic limit, q = 0, Wj becomes aa^ = 1 
and a, are ladder operators. For x= 1> as in |10|, one 
can see using M) that the matrix product (||) is equivalent 
to a problem of counting ID random walks that do not 
return to the origin J5j. For general q, the evaluation of 
(^) poses a q-combinatoric problem, the solution of which 
we now outline. 

We take advantage of recent techniques and results C3] 
for the calculation of matrix products such as (||). The 
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approach is based on the fact that the eigenstates of the 
operator x = a+a^ (analogous to the position operator in 
the usual harmonic oscillator) can be expressed in terms 
of g-deformed Hermite polynomials, whose orthogonality 
properties and generating functions are known. Decom- 
posing (W\ and \V) onto the eigenbasis of x allows an 
integral representation of F(N). From this expression 
the large N behavior can be extracted by using standard 
asymptotic analysis detailed in [Q . The results are sum- 
marized in Table II in which 



K = 



(l-q n ) 4 



(10) 



Using Eq. ^|and the results of Table II, one can calculate 



X < g 
x = q 2 

g 2 <X<l 
X = l 

X > i 
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TABLE II. Explicit expressions of the parameters in the 
generic decay form F(N) - = A(l - B) N+J N" 1 . For 
X < 1, F°° = and far X > 1, F, = l- fll la- 
the asymptotic density ([I]) for any initial spatial distribu- 
tion of particles. In Table I we present the results found 
with G(N; X) given by the Poisson initial condition. 

In summary we have studied stochastic ballistic anni- 
hilation and coalescence in ID. The asymptotic density 
was calculated exactly for Poisson initial conditions. The 
resulting phase diagram is described by two parameters 
q and \- The first is a measure of stochasticity while the 
second encodes information about the reaction processes 
and the initial densities. The phase diagram contains two 
new regimes which were not known before. 

One application of our results is to generalize the sur- 
face growth model of g. In that model down (up) steps 
of a ID interface move dctcrministically to the right (left) 
and annihilate on meeting. Thus the surface smoothcns 
with time. Our results for tjl = f]R = allow one to con- 
sider the effect of a probability q of a new terrace being 
formed when steps meet and the effect of an overall tilted 
surface (x 7^ 1). We see that for % = 1 the probability q 
does not affect the long-time smoothing of the interface, 



whereas for a tilted interface the smoothing behaviour 
changes for q 2 > \- 

Our results are exact in the long-time limit. It would 
be interesting to study the approach to this asymptotic 
behavior. One way to do this would be through extensive 
numerical simulations although, as the cross-over time 
is expected to grow with q, the simulation time needed 
can be very large. Previous numerical studies Jl6| were 
performed at relatively short times and indicated that 
along the line \ = \ the decay exponent depends on q 
which is not the case in the true long-time limit. 

Our results allow extension to other initial spatial dis- 
tributions by using the results in Table II and appropriate 
forms for G(N,X) in (||). We expect other distributions, 
for which the number of particles in a macroscopic re- 
gion obeys a central limit theorem, to exhibit the same 
phases. However, power law distributions might generate 
different behavior. 

Generalizations to higher dimensions and more than 
two velocities are desirable. Indeed, even in one dimen- 
sion and with reaction probability one, models with more 
than two velocities have been shown to exhibit rich be- 
havior fn] ] . It would be interesting to try and generalize 
the analytical approach presented here to that case. 
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